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Abstract. Let K := {x : g(x) < 1} be the compact sub-level set 
^ . of some homogeneous polynomial g. Assume that the only knowledge 

about K is the degree of g as well as the moments of the Lebesgue mea- 
sure on K up to order 2d. Then the vector of coefficients of g is solution 
of a simple linear system whose associated matrix is nonsingular. In 
other words, the moments up to order 2d of the Lebesgue measure on K 
encode all information on the homogeneous polynomial g that defines 
K (in fact, only moments of order d and 2d are needed). 



1. Introduction 

The inverse problem of reconstructing a geometrical object K C W 1 from 
the only knowledge of moments of some measure \x whose support is K is a 
fundamental problem in both applied and pure mathematics with important 
applications in e.g. computer tomography, inverse potentials, signal process- 
ing, and statistics and probability, to cite a few. In computer tomography, 
for instance, the X-ray images of an object can be used to estimate the mo- 
ments of the underlying mass distribution, from which one seeks to recover 
the shape of the object that appears on some given images. In gravimetry 
applications, the measurements of the gravitational field can be converted 
into information concerning the moments, from which one seeks to recover 
the shape of the source of the anomaly. 

Of course, exact reconstruction of objects K C M n is in general impossible 
unless K has very specific properties. For instance, if K is a convex polytope 
then exact recovery of all its vertices has been shown to be possible via a 
variant of what is known as prony method. Only a rough bound on the 
number of vertices is required and relatively few moments suffice for exact 
recovery. For more details the interested reader is referred to the recent 
contribution of Gravin et al. [2] and the references therein. On the other 
hand, Cuyt et al. [I] have shown that approximate recovery of a general 
re-dimensional shape is possible by using an interesting property of multi- 
dimensional Pade approximants, analogous to the Fourier slice theorem for 
the Radon transform. 
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Contribution. From previous contributions and their references, it is trans- 
parent that exact recovery of an n-dimensional shape is a difficult problem 
that can be solved only in a few cases. And so identifying such cases is of 
theoretical and practical interest. The goal of this paper is to identify one 
such case as we show that exact recovery is possible when K C 1™ is the 
(compact) sublevel set {x G R ra : <?(x) < 1} associated with an homoge- 
neous polynomial g. By exact recovery we mean recovery of all coefficients 
of the polynomial g. In fact, exact recovery is not only possible but rather 
straightforward as it suffices to solve a linear system with a nonsingular ma- 
trix! Moreover, only moments of order d and 2d of the Lebesgue measure on 
K are needed. As already mentioned, exact recovery is possible only if K 
has very specific properties and indeed, crucial in the proof is a property of 
levels sets associated with homogeneous polynomials (and in fact, also true 
for level sets of positively homogeneous nonnegative functions). 

2. Main result 

2.1. Notation and definitions. Let R[x] be the ring of polynomials in the 
variables x = (x\,. . . , x n ) and let R[x]<j be the vector space of polynomials 
of degree at most d (whose dimension is s(d) := ( n ,j^))- For every d G N, 
let := {a G N n : \a\ (= = d}, and let v d (x) = (x°), a G N n , be 
the vector of monomials of the canonical basis (x Q ) of Rbx:]^. Denote by St 
the space of k x k real symmetric matrices with scalar product (B,C) = 
trace (BC); also, the notation B >z (resp. B y 0) stands for B is positive 
semidefinite (resp. positive definite). 

A polynomial / G R[x]d is written 

x-»/(x) = 

for some vector of coefficients f = (f a ) G W^ d \ 

A real-valued polynomial g : R n — > R is homogeneous of degree d (d G N) 
if flf(Ax) = \ d g(x) for all A and all x£l. Given g G R[x], denote by G C W 1 
the sublevel set {x : g(x) < 1}. 

If g is homogeneous then G is compact only if g is nonnegative on W n 
(and so d is even). Indeed suppose that g(xo) < for some xo G W 1 ; then 
by homogeneity, g(Axo) < for all A > and so G contains a half-line and 
cannot be compact. 

2.2. Main result. The main result is based on the following result of inde- 
pendent interest valid for positively homogeneous functions (and not only 
homogeneous polynomials). A function / : W 1 — > R is positively homoge- 
neous of degree d G R if /(Ax) = A d /(x) for all A > and all x G R n . 

Lemma 1. Let f : R n — > R be a measurable, positively homogeneous and 
nonnegative function of degree < d G R, with bounded level set {x : /(x) < 
1}. Then for every k G N and a £ N": 
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(2.1) / x"/(x) fc dx = — — — . / x a dx. 

J{x:/(x)<l} n + ™ + M ./{x:/(x)<l} 



Proof. To prove (12. ip we use an argument already used in Morosov and 
Shakirov 13 0]. With a G N n , let d := (a 2 ,...,a n ) G N™^ 1 and define 
z := (z 2 , . . .,z n ). 

Let : M + ->Kbe measurable and consider the integral J K „ (f)(g(x)) x Q dx. 
Using the change of variable x\ = t and ccj = tZi for all i = 2, . . . , n, and 
invoking homogeneity, one obtains: 

/ </>(/(x))x Q <ix = / <f>(t d f(l,z 2 ,...,z n ))t n+ W- 1 z & d(t,z) 
= d~ l ( f°° u^+^y^^du) x A a 



with A* = [ z^/(l,z)"( n+ l Q l)/ d dz. 
Hence the choices t h-> </>(£) := I[o,i] (^) and t := i fc I[ 0i i](i) yield 

d/ X Q dx = A a f 1 V (n+\a\)/d-l du= ^ad 

J{x:/(x)<i} Jo n + [a| 

A n d 



d [ /(x) fc X°dx = A a f u (.n+kd+\a\)/d-l du 

J{x:/(x)<l} JO 



n + kd+ \a\ 



respectively. And so (|2.ip follows. □ 

With g G M[x]d being an homogeneous polynomial of degree d, consider 
now the matrix M<i(A) of moments of order 2d associated with the Lebesgue 
measure on G = {x : <?(x) < 1}, that is, Md(A) is a real square matrix with 
rows and columns indexed by the monomials x a , a G NJJ, and with entries 

(2.2) M d (X)[a, /3] = [ x Q+/3 dx =: \ a+f} , Va,/3eN2. 

■/{x:g(x)<l} 

So all entries of M^(A) are moments of order 2d. Our main result is as 
follows: 

Theorem 2. Let g G M[x]d be homogeneous of degree d with unknown co- 
efficient vector g € and with compact level set G = {x : g(x) < 1}. 
Assume that one knows the moments A = (A a ) for the Lebesgue measure on 
G, for every a G N n with \a\ = 2d and \a\ = d. Then: 

(2-3) g = ^M d (A)- 1 AW 

where A^ = (A a ), a G Nj, is f/ie vector of all moments of order d. 
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Proof. Use (EjJ) with fc = 1 and 



a- 



j3eN n : \p\=2d 

or in matrix form 

M d (A)g 



= d to obtain 
n + d 

An., 



n + 2d 



V lal = d, 



n + d 



A„ 



ra + 2d 

from which the desired result follows if M^(A) is non singular. But this 
follows from the fact that G has nonempty interior. □ 

There are alternative ways for obtaining g from the moments A. It suffices 
to apply (|2.ip for a family T of multi- indices a 6 N n whose cardinal 
matches the dimension ( n+ rf _1 ) of the vector g. 

For instance, with n = 2 and d = 2 (g is a quadratic form with vector of 
coefficients g = (920)511)502)), g can also be obtained by: 





A20 


Aii 


A02 


-1 


n \ 
H+2 A 00 


g = 


A 30 


A21 


A12 




n+1 \ 
^+3 A 10 




A21 


A12 


A 30 




n+1 a 
L 7I+3 A 01 J 



provided that the above inverse matrix exists. 

Similarly, with n = 2 and d = 4 (g is a quartic form with vector of 

coefficients g = (g i0 , g 3 i, 922, 913, 904)), g can also be obtained by 

-1 r n \ 

^+4 A 00 
n+1 \ 



A 40 


A 31 


A 22 


A 13 


A 04 


A 50 


A41 


A32 


A23 


Al4 


A41 


A32 


A23 


Al4 


A 05 


A 60 


A5I 


A42 


A 33 


A24 


A42 


A 33 


A24 


A15 


A 06 



n+1 

n+5 ' 
n+2 
n+6 ' 
n+2 

n+6 7 



01 



V)2 



provided that the above inverse matrix exists. 
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